Hypersensitivity to perturbations of quantum- chaotic wave-packet dynamics 
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We re-examine the problem of the "Loschmidt echo", which measures the sensitivity to per- 
turbation of quantum chaotic dynamics. The overlap squared M(t) of two wave packets evolv- 
ing under slightly different Hamiltonians is shown to have the double-exponential initial decay 
tx exp(— constant x e 2A °*) in the main part of phase space. The coefficient Ao is the self- averaging 
Lyapunov exponent. The average decay M tx e _Ali is single exponential with a different coefficient 
Ai. The volume of phase space that contributes to M vanishes in the classical limit Ti — > for times 
less than the Ehrenfest time te = ^A^ 1 1 In h\. It is only after the Ehrenfest time that the average 
decay is representative for a typical initial condition. 



PACS numbers: 05.45Mt, 03.65.Sq, 03.65.Yz, 05.45.Pq 



Chaos in classical mechanics is characterised by an ex- 
ponential sensitivity to initial conditions: The separation 
of two trajectories that are initially close together in- 
creases in time cx e xt with a rate given by the Lyapunov 
exponent A. There is no such sensitivity in quantum me- 
chanics, because the overlap of two wave functions is time 
independent. This elementary observation is at the ori- 
gin of a large literature (reviewed in a textbook PJ) on 
quantum characterisations of chaotic dynamics. 

One particularly fruitful line of research goes back to 
the proposal of Schack and Caves Q , motivated by earlier 
work of Peres || , to characterise chaos by the sensitivity 
to perturbations. Indeed, if one and the same state tpo 
evolves under the action of two different Hamiltonians H 
and H + SH, then the overlap 



M(t) = \{^\e %i - H+m ^/ h e^ Ht / h \^ 



(1) 



is not constrained by unitarity. Jalabert and Pastawski[|j 
discovered that M(t) (which they referred to as the 
"Loschmidt echo") decays cx e~ At if ipo is a narrow wave 
packet in a chaotic region of phase space, providing an ap- 
pealing connection between classical and quantum chaos. 

The discovery of Jalabert and Pastawski gave a new 
impetus || to what Schack and Caves called "hypersen- 
sitivity to perturbations" of quantum chaotic dynamics. 
The present paper differs from this body of literature in 
that we consider the statistics of M(t) as tpo varies over 
the chaotic phase space. We find that the average de- 



cay M(t) oc 



is due to regions of phase space that 



become vanishingly small in the classical limit fi c s — * 0. 
(The effective Planck constant h e g = Ti/Sq is set by the 
inverse of a typical action Sq.) The dominant decay is a 
double exponential oc exp(— constant x e 2At ), so it is truly 
"hypersensitive" . The slower single exponential decay is 
recovered at the Ehrenfest time te = A 1 1 ln7i e ff|. 

Before presenting our analytical theory we show in 
Fig. 1 the data from a numerical simulation that illus- 
trates the hypersensitivity mentioned above. The Hamil- 



tonian is the quantum kicked rotator jjj 
TI : : — +K cosx^5(t-5 



?! 
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^ i dx 



(2) 



The perturbed Hamiltonian H' = H + SH is obtained 
by the replacement K — > K + 5K . The coordinate x is 
periodic, x = x + 2n. To work with a finite dimensional 
Hilbert space we discretize x k = 2irk/N, k = 1,2,. ..TV. 
The momentum p m — mfi is a multiple of h, to ensure 
single-valued wavefunctions. For h = h c g = 2tt/N the 
restriction to the first Brillouin zone results in a single 
band p m = 2i:m/N, m — 1,2, ...N. The time evolution 
e -iHn/h = jjn a fter n periods, of the initial Gaussian 
wave packet ipk = N~ x / 2 exp {ttN~ x [2im,ok — (k — fco) 2 ]) > 
is given by the Floquet operator in the ^-representation: 



1 finCk'-k) 2 NK 2irk\ 
Uk , k = —e W ^eos— j 



(3) 



We use the fast Fourier transform algorithm to compute 
U n for TV up to 10 6 §. 

We study the statistics of M(t) by comparing in Fig. 1 
three different ways of averaging over initial positions 
(mo, fco) of the Gaussian wave packet. We used K = 10, 
SK = 1.6 • 10_^ and N = 10 6 (^ off = 6.28 ■ 1CT 6 ). While 
the average M decays exponentially, the two logarith- 
mic averages have a much more rapid initial decay. We 
estimate that M < 1CT 23 at n = 3 for about 30% of 
randomly chosen initial conditions. For the same point 
n = 3 only 9% of initial conditions (corresponding to 
M > 0.2) account for 80% of the total value of M. The 
typical decay of M(t) is therefore much more rapid than 
the exponential decay of the average M. 

Statistical fluctuations also affect the decay rate 
of M, set by the Lyapunov exponent according to 
Ref. H . The definition of the Lyapunov exponent Ao = 
limt^oo t^ 1 In \6x(t)/Sx(0)\ gives Ao = 1.65 for the classi- 
cal kicked rotator with K = 10. However, M(t) in Fig. 1 
has exponent Ai = 1.1, defined by 



A, 



lim (jt)- 1 ln|fe(t)/()a;(0)|-J. 



(4) 



2 




10 12 14 16 



FIG. 1: The overlap M at t — n for the quantum kicked ro- 
tator for three different ways of averaging: o M, o exp(lnM), 
• exp[- exp(ln(-lnM))]. We took K = 10, SK = 1.6 ■ 1(T 3 , 
N — 10 6 . Averages are taken over 2000 random initial condi- 
tions of a Gaussian wave packet. The dotted line shows the 
Lyapunov decay oc e _nAl with Ai = 1.1. At n = 3 we have 
only an upper bound for the logarithmic averages because 
cancelations in the calculation limit the accuracy. 



Since fluctuations of t^ 1 In \Sx(t)/Sx(0)\ decrease like 
t~ x l 2 the Lyapunov exponent Ao is self-averaging [Q, 
while the Aj's are not. 

For an analytical description we start from the Gaus- 
sian one-dimensional wave packet 



^(• T ) = (^) 1 ex P + ttP 



The wave packet is centered at the point xo(t),p (t) 
which moves along a classical trajectory. Initially (3(t = 
0) = and a(t = 0) = 1. Divergence of trajectories 
leads to the exponential broadening of the packet, thus 
a(t) oc cxp(— 2At). Since a -C 1 for t ^> 1/A, the wave 
packet in ph ase space be comes highly elongated with 
length Z|| = y/h(l + (3 2 )/a and width l± — h/lu. The pa- 
rameter (3 — Ap/ Ax represents the tilt angle of the elon- 
gated wave packet The Gaussian approximation (0) 
breaks down at the Ehrenfest time te = \\~ l \hxW\, 
when in becomes of the order of the size of the system. 

We assume that ip evolves according to Hamiltonian 
H(K) and ip' according to H' = H(K+5K). The overlap 
M = | (ip'\if>) | 2 of the two Gaussian wave packets is 



M = 



45(3- 



■ exp 



a(Sp — (38x) 2 aa'Sx' 
2(a 2 +4S(3 2 )h mh 



(6) 

in terms of the (weighted) mean a = (a + oc')/2, 
(3 = (f3' a + (3a') /(a + a') and difference Sp = p' — pa, 
5x = x' — xq, Sa = a' — a, 5/3 = (3' — (3. In order of 
magnitude, 5j3/ f3 ~ 8a /a ~ 8K -c 1. The displacement 
vector (8x,Sp) has component Aii ~ 5Ke xt parallel to 
the elongated wave packets and component A^ ~ 5K 
perpendicular to them (see Fig. 2). 



Depending on the strength of perturbation one may 
distinguish three main regimes: SK <?i,h< SK < Vh, 
and 5K > y/h. We will consider in detail the interme- 
diate regime Ti < SK < \f% and discuss the two other 
regimes more briefly at the end of the paper. (The simu- 
lations of Fig. 1 are at the upper end of the intermediate 
regime, since SK = 1.6 • 10~ 3 and \fh = 2.5 • 10 -3 .) 
The three regimes may be characterised by the relative 
magnitude of the Ehrenfest time te and the perturbation 
dependent time scale To = i A - 1 1 ln<5i^|. In the interme- 
diate regime one has ^te < tq < te- 

To estimate the relative magnitude of the two terms in 
the exponent of Eq. (||) we write 

5p-P5x= (l+^V^Ax = fSK, (7) 

7T- „2Ai 



a 2 + 4,5(3- 



= Q = 



1 + (ge 2Xt 5K) 2 ' 



(8) 



Here / and g are functions of order unity of time t and 
the initial location Xi,pi of the wave packet. The second 
term in the exponent @ is of order aSx 2 /h ~ SK 2 /%, 
while the first term is of order Q5K 2 /h. Since Q 3> 1 for 
t < 2tq, and since 2r > te in the intermediate regime, 
we may neglect the second term relative to the first term 
within the entire range t < te of validity of the Gaussian 
approximation. Eq. m) thus simplifies to 



M = (SQ) 1/2 eyip[-^Q(fSK) 2 /K 



1 



(9) 



We seek the statistics of M(t) generated by varying 
Xi,pi. The statistics is non-trivial because fluctuations 
in / of order unity cause exponentially large fluctuations 
in M if Q5K 2 jfi 1, which is the case for 2r — te < 
t < te- The average of M is then dominated by the 
nodal lines x n (p) in phase space at which / vanishes (at 
a particular time t). If Ax n is the typical spacing of these 
lines at constant p, then the derivative df jdxi at x n is 
of order 1/Ax n . This yields 



M 



= (aQ) 1/2 



dx 



exp 



A3 



QSK 2 
2h 



{afi/SK 2 ) 1 / 2 = (Vh/SK)e- xt . 



(10) 



Assuming independent fluctuations in the (perturbation 
dependent) distribution of nodal lines and in the rate 
of divergence of trajectories for the individual Hamilto- 
nian, we incorporate fluctuations in A in Eq. (^) via 
cxp(— Xt) — > exp(— Aii), in accordance with ([I]). Hence 
we recover the exponential decay of the Loschmidt echo 
, although with the exponent Ai instead of Ao (in agree- 
ment with the numerics of Fig. 1). The exponential decay 
sets in for t > 2tq — te, while for shorter times M remains 
close to unity . 

The volume V of phase space near the nodal lines con- 
tributing to M is of order V = (Ti/QSK 2 ) 1 / 2 . This vol- 
ume decreases exponentially in time for t < To, reaching 
the minimal value Vq = \JXijSK <C 1 at tq. For larger 




FIG. 2: Schematic illustration of two perturbed wave packets 
in phase space for t < to. 



times V increases, saturating at a value of order unity 
at te ■ We therefore conclude that the average M is only 
representative for the typical decay if t > te- For smaller 
times the average is dominated by rare fluctuations that 
represent only a small fraction of the chaotic phase space. 

To obtain an average quantity that is representative 
for a typical point in phase space we take logarithmic 
averages of Eq. (0) . For t < tq one has 



InM ~ ~(SK 2 /h) exp(2A_ 2 i), 



(11) 



lnln(l/M) = 2X t - hi{h/5K 2 ) + 0(1). (12) 



(The coefficient A_ 2 m InM appears because we aver- 
age the square of displacement.) The double logarithmic 
average (|l2|), given by the self-averaging Lyapunov ex- 
ponent Ao, is least sensitive to fluctuations and is repre- 
sentative for the main part of phase space. The typical 
overlap thus has the double-exponential decay 

M ~ exp ( — constant x (5K 2 /h)e 2Xot ) , (13) 

down to a minimal value Mq ~ exp(— 5K/K) at t = tq. 

The initial decay (13) for t <C to is the same as ob- 
tained in Ref. |l(J for the classical fidelity (defined as the 
overlap of two classical phase space densities). In that 
problem the role of h is played by the initially occupied 
volume of phase space. A superexponential decay of the 
classical fidelity has also been obtained by Eckhardt jll] . 

The origin of the decay ( |l3| ) is illustrated in Fig. 2. For 
t < tq the wave packets are nearly parallel (<5/3 <C a) , dis- 
placed laterally by an amount Aj_ cx SK . Their overlap 
is an exponential function oc exp(— Aj_/ZjJ, where the 
width Zj_ of each wave packet decreases exponentially in 
time cx e~ At . Hence the double-exponential decay. 

For t > To (when 5(3 ^> a) the overlap of the two wave 
packets is dominated by their crossing point x c ,p c . The 
overlap M ~ exp (— constant x \x c — Xq\ /I 2 ) now increases 

with time, because Zi 



. Since \x c — xq\ ~ A±/5(3 ■ 



/, the crossing point falls outside of the range of validity 
of the Gaussian approximation unless |/| <C 1. The result 
( |To| ) is justified (because it is dominated by nodes of /), 
but we can not use the Gaussian approximation to extend 
the formula ([ll]) for the typical decay to t > tq . The typ- 
ical decay and the average decay become the same at te, 




FIG. 3: Two perturbed wave packets in phase space for 
te < t < te + 2t . The lines show p(x) (solid) and p'(x) 
(dashed) extracted from the Husimi function evolved with 
the quantum kicked rotator, for TV" = 10 6 , K = 7, SK = 0.1, 
n — 5. Dots show the crossing points Xj that contribute to 
the overlap in stationary phase approximation. 



so the typical M should increase from its minimal value 
M at r to the value M E = {Vh/8K)e~ XTE = h/SK at 
te- Both Mo and Me are <C 1, but Mo is exponentially 
small in 5K/h while Me is only algebraically small. 

For t > te one can use the semiclassical WKB de- 
scription of elongated wave packets, along the lines of 
Refs. jl^]. The phase space representation of the wave 
function ip is concentrated along the line on the torus 
p(x) of length l\\ ~ Vhe xt ^> 1, see Fig. 3. The function 
p(x) is multivalued and each branch k has a WKB wave 
function with amplitude pk w and phase o^: 



i<y k /h 



Pk = da k /dx. 



(14) 



For SK > Ti the overlap of two oscillating wave functions 
ip,tp' of the form ([ll]) may be found in stationary phase 
approximation. The stationary points are given by the 
crossings p{xj) — p'(xj) of the two lines p(x), p'{x) given 
by the evolution with Hamiltonians H, H' . For te < t < 
te + 2to the number of crossings N c is proportional to Zy 
and independent of SK. This is because both the lateral 
displacement of p and p' and their relative angle are of 
the same order SK. (In Fig. 3 we have Zy ~ 20iV c .) Each 
crossing contributes to (f/'IV'') an amount 



P, 



p(xj)p'(xj) / cZxexp 



+ ic 



(15) 



where k = d 2 (a — a')/dx 2 \ x . ~ SK and h<pj = cr(xj) — 
a'(xj). The phase <pj varies randomly from one crossing 
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FIG. 4: Decay of the average overlap for the quantum kicked 
rotator (K = 10) in the golden-rule regime. We keep T = 
0.023(/V5#)^_fixed by taking SK = 1/TV. Circles from bottom 
to top give M for TV = 10 3 , 10 4 , 10 5 , 10 6 ,4 • 10 6 . The inset 
demonstrates that no scales with log TV, as expected for the 
Ehrenfest time. 

to the other, leading to 

M= h V e '(^-4y) 5 —e~ Xlt (16) 

l\8K V ~ U5K- 5K e ■ [W) 

II j,f=i " 

Because of the large number of crossings there is now 
little difference between M and logarithmic averages. For 
t > te + 2to the number of crossings becomes N c ~ SKI 2 . 
(The distance between almost parallel segments of p'(x) 
is of order and the line p(x) crosses at the angle SK 
about SKlu 3> 1 segments per unit length.) This leads 
to saturation of the overlap at M ~ h. 

This completes our discussion of the intermediate 
regime Ti < SK < yh. We conclude with a brief dis- 
cussion of the two other regimes. For SK > VK the 
longitudinal displacement of the packets exceeds their 
lengths A || > The logarithmic averages now remain 



the same, but M is changed. The dominant contributions 
to M are now given by the rare events for which both 
A_l and A|| vanish. This leads to M ~ (h/SK 2 )e~ Xlt 
for t < 2tq. (The same Lyapunov decay as in the in- 
termediate regime, but with a much smaller prefactor.) 
For 2tq < t < te the length of each wave packet remains 
small, Z|| ~ y/he xt <C 1, but the displacement saturates 
at the maximal value A|| ~ 1. In this time range the 
average overlap has a plateau at M ~ U/SK. Finally, for 
t > t e the decay © M ~ {^JK/8K)e' Xlt is recovered. 

In the remaining regime SK <C h we find from Eq. @ 
that M(t) remains close to unity for t < te, regardless of 
the initial location of the wave packet. (This also results 
in insensitivity to the way of averaging.) The golden-rule 
decay with rate T ~ (5K/K) 2 , sets in only after the 
Ehrenfest time: M ~ exp[— T(t — te)] for t> t E - These 
results are depicted in Fig. 4. Golden-rule decay persists 
until the Heisenberg time tn — l/h or the saturation 
time r _1 |ln?i|, whichever is smaller. (Only the initial 
decay is shown in Fig. 4.) The Gaussian decay H sets in 
for t > t H provided that SK < h 3/2 . 

In summary, we have shown that statistical fluctu- 
ations play a dominant role in the problem of the 
Loschmidt echo on time scales below the Ehrenfest time. 
While the decay of the squared overlap M(t) of two per- 
turbed wave packets is exponential on average, as ob- 
tained previously |Q, the typical decay is double expo- 
nential. It is only after the Ehrenfest time that the main 
part of phase space follows the single-exponential de- 
cay of M. The Ehrenfest time has been heavily studied 
in connection with the quantum-to-classical correspon- 
dence |^| . The role which this time scale plays in sup- 
pressing statistical fluctuations has not been anticipated 
in this large body of literature. 
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